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We discuss new techniques developed in recent years for performing one-loop calculations in QCD, and present 
an example of results from the process — > Vqqgg. 



1. INTRODUCTION 

Perturbative QCD, and jet physics in particu- 
lar, have matured sufficiently that rather than be- 
ing merely subjects of experimental studies, they 
are now tools in the search for new physics. This 
role can be seen in the search for the top quark, 
as well as in recent speculations about the im- 
plications of supposed high- .Ex deviations of the 
inclusive-jet differential cross section at the Teva- 
tron. One of the important challenges to both 
theorists and experimenters in coming years will 
be to hone jet physics as a tool in the quest for 
physics underlying the standard model. As such, 
it will be important to measurements of parame- 
ters of the theory or of non-perturbative quanti- 
ties such as the parton distribution functions, as 
well as to searches for new physics at the LHC. 

Jet production, or jet production in associa- 
tion with identified photons or electroweak vector 
bosons, appears likely to provide the best infor- 
mation on the gluon distribution in the proton, 
and may also provide useful information on the 
strong coupling a s . 

In order to make use of these final states, we 
need a wider variety of higher-order calculations 



of matrix elements. Indeed, as we shall review 
in the next section, next-to-leading order calcu- 
lations are in a certain sense the minimal useful 
ones. On the other hand, these calculations are 
quite difficult with conventional Feynman rules. 
In following sections, we will discuss some of the 
techniques developed in recent years to simplify 
such calculations. 



NEXT-TO-LEADING 
CULATIONS 



ORDER CAL- 
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While leading-order calculations often repro- 
duce the shapes of distributions well, they suf- 
fer from several practical and conceptual prob- 
lems whose resolution requires the use of next-to- 
leading order calculations. These problems are 
tied to the various logarithms that can arise in 
perturbation theory. 

The first of these logarithms are 'UV ones, con- 
nected with the renormalization scale. We are 
forced to introduce a renormalization scale \x in 
order to define the coupling, a s (/Lt); but physi- 
cal quantities, such as cross sections or differ- 
ential cross sections, should be independent of 
/Lt. When we compute such a quantity in pertur- 
bation theory, however, we necessarily truncate 
its expansion in a s , and this introduces a spu- 
rious dependence on /j,. This dependence is sig- 
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nificant in real-world applications of perturbative 
QCD because the coupling is not that small, be- 
cause it runs relatively quickly, and because we 
are interested in processes with a relatively large 
number of colored 'final'-state partons. Together, 
these effects can lead to anywhere from a 30% 
to a factor of 2-3 normalization uncertainty in 
predictions of experimentally-measured distribu- 
tions. At leading order, the only dependence on fx 
comes from the resummation of logarithms in the 
running coupling a s (fi), and the scale choice is 
arbitrary. At next-to-leading order, however, the 
virtual corrections to the matrix element intro- 
duce another dependence on /i. This dependence 
can — and in practice, often does — reduce the 
over-all sensitivity of a prediction to variations in 
/i. (I should stress, however, that while varying fi 
by some preset amount, say a factor of two up and 
down from a typical scale, gives an indication of 
the sensitivity of or uncertainty in the calculation, 
it does not give an estimate of the error involved; 
for that one needs a next-to-next-to-leading order 
calculation.) 

The other logarithms are the 'IR' ones, con- 
nected with the presence of soft and collinear 
radiation. Jets in a detector are not infinitely 
narrow pencil beams; they consist of a spray of 
hadrons spread over a finite segment. Experi- 
mental measurements of jet distributions and the 
like depend on resolution parameters, such as the 
jet cone size and minimum transverse energy. In 
a leading-order calculation, jets are modelled by 
lone partons. As a result, these predictions don't 
depend on these parameters (or have an incorrect 
dependence on them). In addition, the internal 
structure of a jet cannot be predicted at all. 

At next-to-leading order, one necessarily in- 
cludes contributions with additional real radia- 
tion, which either shows up inside one of the jets, 
or as soft radiation in the event. This introduces, 
for example, the required logarithmic dependence 
on the cone size AR. 

At a more conceptual level, we must remember 
that jet differential cross sections are multi-scale 
quantities, involving not only a hard scale (say a 
jet transverse energy Et) but also scales charac- 
terizing the resolution defining a jet (for exam- 
ple, EtAR or -E-rmin)- As a result, the pertur- 



bative expansion is not one in a s alone, but con- 
tains logarithms and logarithms-squared of ratios 
of scales. These logarithms, which arise from the 
infrared structure of gauge theories, might spoil 
the applicability of perturbation theory if they 
grow too large. Only a next-to-leading order cal- 
culation can tell us if we are safe. 



3. ORGANIZING THE CALCULATION 
OF AMPLITUDES 

A leading-order calculation of an n-jet pro- 
cess at a hadron collider such as the Tevatron 
or LHC requires knowledge of the parton distri- 
bution functions of the proton; of a s ; and of the 
tree-level matrix element for the 2 — > n parton 
process. A next-to-leading order calculation of 
the same process requires in addition the tree- 
level matrix element for the 2 — > n + 1 parton 
process, and the one-loop matrix element for the 
2 — * n parton process. (In addition, it requires 
a general formalism such as that of refs. Jl|-Ql for 
cancelling the infrared singularities analytically 
while allowing a numerical calculation of fully- 
differential observables.) 

It is the calculation of the one-loop matrix ele- 
ment which presents the greatest difficulty in the 
traditional Feynman-diagram approach. This dif- 
ficulty arises from the enormous number of dia- 
grams, the large amount of vertex algebra in each 
diagram, and the complexity of loop integrals 
with many powers of the loop momentum in the 
numerator. A brute-force approach might easily 
lead to expressions thousands of times larger than 
an appropriate representation of a result. 

The first step in developing a more efficient 
method for performing these calculations is to 
take advantage of the lessons from earlier develop- 
ments in tree-level calculations, and from string- 
based methods for one-loop calculations. These 
include 

• Color decomposition: one should write the 
amplitude as a sum over color factors mul- 
tiplied by color-ordered subamplitudes, and 
calculate these latter coefficients. For a one- 
loop amplitude for — > Vqqg ■ ■ ■ g, the color 
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decomposition takes the form 

i=i <Tes„_ 2 /s„- 2 ;i 
x^(l„2^(3...n);y) (1) 

where Gr^ (3 • • • n) is the leading-color 
trace, equal to the number of col- 
ors N c times the tree-level color factor 
(T a3 ■ ■ ■ T an Y i 1 2 . The remaining color factors 
are subleading-color ones, 

Gr™- >1 (3---n) = Tr(T° 3 • • • T^ 1 ) 

x (t%'+2 . ,.r")| . (2) 

• Spinor helicity method Q: one should cal- 
culate helicity amplitudes (and square or 
compute interference terms only numeri- 
cally), using spinor products (ij) and [ij]. 

• Use of supersymmetric decompositions and 
supersymmetry identities, where applica- 
ble. 

• Use of permutation identities j|] that allow 
rewriting subleading-color amplitudes as a 
sum of permutation of certain leading-color 
'primitive' amplitudes. One should then 
calculate only primitive amplitudes, which 
correspond to leading-color amplitudes with 
a definite clockwise or counter-clockwise 
orientation of internal fermion lines, and in 
which the external fermion legs are not nec- 
essarily adjacent. 

The reader will find a review of the techniques 
described below in ref . . 

4. PRIMITIVE AMPLITUDES 

In principle, one may imagine calculating prim- 
itive amplitudes using color-ordered Feynman 
rules. One can do better, calculating them us- 
ing a combination of string-based rules and rules 
inspired by string theory. Can one do better yet? 
The unitarity-based technique, combined with ju- 
dicious use of the collinear limits, provides an 
even better way of calculating these quantities. 



In general, a one-loop amplitude has absorptive 
pieces; the basic idea behind the unitarity-based 
technique is that the cut-containing terms in the 
amplitude can be determined by dispersion rela- 
tions from the dispersive parts, which are in turn 
related to a product of tree amplitudes. As I will 
review below, in supersymmetric theories, the ba- 
sic technique in fact determines cut-free pieces as 
well. 

In calculating loop diagrams, we use dimen- 
sional regularization for both the ultraviolet and 
infrared singularities; at one loop, it is possible 
to choose a variant (dimension reduction) that 
preserves supersymmetry manifestly. (The trans- 
lation to the conventional scheme is detailed in a 
paper by Kunszt, Signer, and Trocsanyi ||.) 

5. INTEGRALS 

Let us consider the calculation of an n-point 
one-loop amplitude in a gauge theory. The calcu- 
lation necessarily involves up to n-point one-loop 
integrals; and given the powers of momentum in 
the three-gluon vertex, may involve up to m pow- 
ers of the loop momentum in any m-point inte- 
gral. Using either a combination of Passarino- 
Veltman || and van Neerven-Vermaseren |]lof re- 
duction techniques, or equivalent ones we 
can rewrite any such integral (in D — 4 — 2e) 
as a linear combination of scalar box integrals, 
scalar three-mass triangles, scalar or tensor one- 
or two-mass triangles, and scalar or tensor bubble 
integrals. If we start with an amplitude where all 
lines are massless, the internal lines of all these 
integrals will also be massless; but the external 
legs will be sums of the original external mo- 
menta, and thus possibly massive. The scalar 
boxes we must consider in general will thus have 
anywhere from one to four external masses, for 
example. We will denote the set of integrals that 
appears in the computation of a n-point ampli- 
tude by J- n . The amplitude can be written as a 
sum over integrals in this set, with coefficients Ci 
which are rational functions of the spinor prod- 
ucts and Lorentz invariants. 
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6. UNIQUE RECONSTRUCTION OF 
AMPLITUDES 

One of the basic tools in the method is the re- 
construction theorem, which tells us under which 
conditions the knowledge of the cuts (in four di- 
mensions) suffices to reconstruct the amplitude 
uniquely. It states that if, for every loop integral 

f d D k Poly(^) ( . 

J (2n)D F(fc- Pl ) 2 ---(fc-Pn-i) 2 { ' 

which appears in the amplitude, the degree of the 
numerator polynomial is less than or equal to n — 
2 (less than or equal to 1 if n = 2), then the 
amplitude is determined uniquely by its cuts in 
D = 4. This does not imply that the amplitude 
is necessarily free of rational terms, but merely 
that in this case, the rational terms are always 
linked to the terms containing cuts. 

Using the background-field method, and exam- 
ining the effective action, one can show that the 
theorem applies to the complete amplitude in su- 
persymmetric theories. In TV = 4 supersymmet- 
ric gauge theories, one can go further, and show 
that in the background-field method, the highest 
power of the loop momentum that can appear in 
the numerator of an n-point integral is n — 4; this 
in turn implies that in these theories, the ampli- 
tude can be expressed as a sum of scalar box in- 
tegrals with rational coefficients a (triangles and 
bubbles don't enter). 

7. SEWING AMPLITUDES 

The key point in the unitarity-based method is 
that we sew tree amplitudes, not tree diagrams. 
This allows us to take advantage of all the cancel- 
lations and reductions in numbers of terms that 
have already occurred in the process of computing 
the tree amplitudes. 

A calculation using the unitarity-based method 
proceeds as follows. We want to compute the co- 
efficients a of the integrals in the set T n - To 
this end, we consider in turn cuts in all possible 
channels. For a given channel, we form the cut 
in that channel, summing over all intermediate 
states; this gives rise to a phase-space integral of 



the form 

d D LiPS(-e 1 ,e 2 ) 

xAt ee (-e 1 ,...,e 2 )A% ee (-e 2 ,...,e 1 ), (4) 

where arc the four-dimensional on-shcll mo- 
menta crossing the cut, and where the A* 1 ™^ are 
the color-ordered tree subamplitudes on the two 
sides of the cut. Using the Cutkosky rules, we 
can rewrite this expression as the absorptive or 
imaginary part of a loop amplitude, 

u 4 tree 1 \ tree 1 /r\ 

wy 5 L * (5) 

Using the power-counting theorem in those cases 
where it applies (or simply up to a rational ambi- 
guity to be fixed later in cases where it does not 
apply), we can recover the real parts by drop- 
ping the subscript "cut". We then perform the 
usual reductions on the resulting loop integral, 
and extract the coefficient of any function in T n 
containing a cut in the given channel. (Functions 
not containing a cut in the given channel should 
be dropped.) Finally, we reassemble the final an- 
swer by considering all channels. 

Certain integrals, such as the box integrals, 
have cuts in more than one channel. Consider- 
ing both channels provides us with a cross check 
the coefficients as computed in both channels 
must agree — or alternatively we could reduce 
the amount of work we must do by considering 
only one channel. The latter choice is particu- 
larly appropriate when computing amplitudes in 
aniV = 4 supersymmetric gauge theory, since as 
noted above all amplitudes in this theory can be 
written in terms of scalar boxes. 

8. RATIONAL TERMS 

In nonsupersymmetric theories (such as QCD), 
there is a remaining rational part which cannot 
be determined in this manner. There are two ap- 
proaches to determining these pieces: the use of 
collinear limits, and extending the cuts to 0(e). 

The collinear approach is based on the uni- 
versal factorization of amplitudes in the collinear 
limit, 

A^(...,a,b)^ 
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A s =± 



+Ci- 1 ° op ( a A ",6^K- c 1 (---;S AE ) 



(0) 



where the splitting amplitudes C are universal 
functions depending only on the collinear mo- 
menta and their helicities A, and on the momen- 
tum fraction z (k^, — k a + fcf,, k a — zk-%). This 
limit is useful for deducing the rational pieces be- 
cause the cut terms alone will not produce the de- 
sired collinear limit; one must add rational terms 
in order to obtain the correct limit. However, 
there is no proof that the terms deduced by this 
method are unique (though it is likely true for 
more than five external legs; for five external legs, 
there is an ambiguity arising from the existence 
of a term which contains no cuts but is collinear- 
finite). Such results thus need to be checked 
(for example, numerically) against results from 
another method. 

Another possible method again uses the cuts, 
but at 0(e) rather than just to 0(e°). The ba- 
sic point is that amplitudes in a massless gauge 
theory have an over-all power of (— s)~ e , where 
s is an invariant. The "rational" pieces there- 
fore do contain cuts at 0(e), and can be de- 
duced by sewing tree amplitudes, where the mo- 
menta crossing the cuts are taken to be on-shell in 
(4 — 2e) dimensions rather than four dimensions. 
For scalars, this is equivalent to computing mas- 
sive rather than massless amplitudes, followed by 
an appropriately weighted integration over the 
"mass" (really the (— 2e)-dimensional component 
of the momentum). 

9. AN ALL-MULTIPLICITY RESULT 

We had previously obtained [fT2f, with Dave 
Dunbar, an explicit formula for the one-loop am- 
plitude, in the N = 4 supersymmetric gauge the- 
ory, with an arbitrary number of external glu- 
ons with the helicity configuration of the Parke- 
Taylor tree- level amplitudes |l3|,[l4|], namely all 
but two gluons carrying the same helicity. (Let 
us take the majority to have positive helicity, 
the two opposite-helicity gluons to have nega- 
tive helicity; exchanging these two amounts to a 
complex conjugation.) Here, we will outline the 



derivation of this result using the unitarity-based 
method. Such results demonstrate the power of 
the method, because they would require the com- 
putation of an infinite number of Feynman dia- 
grams in the traditional approach. 

This amplitude is uniquely determined by its 
cuts, and can be expressed entirely in terms of 
scalar box integrals. There are two kinds of cuts 
to consider, one where both opposite-helicity glu- 
ons are one the same side of the cut, the other 
where one such gluon appears on each side of 
the cut. The former cut receives contributions 
only from gluons crossing the cut, the latter from 
fcrmions and scalars as well; but it turns out that 
after summing over the N = 4 supermultiplet 
(using supersymmetry Ward identities p5[), the 
structure of the second type of cut is similar to 
that of the first. It thus suffices to consider the 
first kind. 

Let us suppose that the opposite-helicity glu- 
ons are on the left side of the cut (drawn vertically 
through the diagram). Then all the external legs 
on the right-hand side of the cut have positive 
helicity; as a result, the tree amplitude on the 
right-hand side will vanish unless both legs cross- 
ing the cut on the right have negative helicity. 
Since we adopt the convention that all momenta 
are directed inwards into an amplitude, the cut- 
crossing legs flip sign as we cross the cut to the 
left. We must flip their helicities as well, which 
tells us that all legs of the tree amplitude to the 
left of the cut have positive helicity, except the 
two external legs carrying negative helicity. 

Thus both tree amplitudes on either side of the 
cut have the helicity configuration of the Parke- 
Taylor amplitudes 13 14J] , so that simple all-n for- 
mulae are available for them. Most of the factors 
in these amplitudes are independent of the mo- 
menta crossing the cut, and so can be pulled out 
in front of the cut integral (Q). Indeed, most of 
these factors reassemble into the tree amplitude 
itself, so that (up to a phase) we obtain 



AT" J d D lAF§ 

((mi - 1) mi) (gi l 2 ) 2 (m 2 (ma + 1)) 



(7) 



N ((mi - 1) lj) (h mi) (m 2 £ 2 ) (li (m 2 + 1)) 
where (mi — 1, mi) and (m 2 , m 2 -I- 1) are the pairs 
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of legs adjacent to the cut. The four factors in the 
denominator give rise to four propagators, which 
along with the two propagators crossing the cut 
means that for an arbitrary number of external 
gluons, we need consider only a hexagon integral. 
Indeed, it turns out to be a rather special hexagon 
integral; rearranging the spinor products in its 
numerator, it decomposes into a sum of one- and 
two-mass scalar boxes. From this decomposition 
and the uniqueness theorem, we then obtain the 
complete one-loop N — 4 supersymmetric ampli- 
tude. 

10. Z DECAY TO FOUR JETS 

The decay of the Z into jets has been stud- 
ied extensively at LEP and SLC. The three-jet to 
two-jet ratio is one of the better methods of de- 
termining a s (Mz), but at present this extraction 
of the running coupling is dominated by theoret- 
ical uncertainties. The resolution of these uncer- 
tainties requires a next-to-next-to-leading order 
(NNLO) calculation of three-jet production. 

Such a calculation will contain as ingredients 
the two-loop corrections to Z — > qqg, which 
have yet to be computed; the Z — > qqggg and 
Z — > qqq'q'g matrix elements at tree level, which 
are known pq; and the one- loop corrections to 
Z — > qqgg and Z — > qqq'q', which we are cur- 
rently computing. (In addition, it will require an 
extension of one-loop techniques for cancellation 
of IR singularities, such as the slicing method of 
Giele and Glover g, to NNLO.) 

These one-loop corrections are also of interest 
outside of this NNLO context, as they can be used 
to produce a program calculating Z — > 4 jets 
at next-to-leading order. This process is the 
lowest-order one in which the gluon and quark 
color charges can be measured independently. A 
next-to-leading order calculation will allow reli- 
able limits to be set, using LEP and SLC data, 
on other possible light colored fermions (such as 
light gluinos). 

We present here the result for one of the 
leading-color subamplitudes of the Z — > qqgg pro- 
cess. It is convenient to tack on the amplitude for 
the Z production from two leptons, so that all 
external legs are massless; in the formulae below, 



these legs will be k$ and k&. Furthermore, we will 
also include a factor of 1/ s^q (corresponding to a 
photon propagator replacing the Z). 

To express this result, we should first record 
the corresponding tree amplitude, 

^4+y(lg )2 + ,3 ,4^;5^,6^) = 

[12] (13) (45) (6+|(l + 2)|3+) 

(1 2) S 23 ii 23 S 56 

(3 4) [2 4] [16] (5- 1 (3 + 4) |2~) 

[34] S23^234S56 

(5-|(3 + 4)|2-) (6+|(l + 2)|3+) 
(1 2)[3 4] S23S56 



(8) 



define several functions (roughly corresponding to 
the various box and triangle integrals that ap- 
pear) 
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ln(r) 



Li(r) = 



Lo(r) + 1 



1 — r 1 — r 

Ls_i(ri,r 2 ) = Li 2 (l - r x ) + Li 2 (l - r 2 ) 

7T 2 
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as well as a flip symmetry, 
flip = {1 <-> 4, 2 <-> 3, 5 <-> 6; ( ) <-> 
The one-loop amplitude is then 



]}■ 



A 1— loop 

A i+V;l ~ 

where 



F n 



(10) 
(11) 



V a - S 



-Sl2 



S23 



-S34 



3 

2~t 



^56 



(12) 
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F g = Bl LS_! 



and where 

S23 



-S34 



"^234 —^234 / 



+ B 5 Ls-l| 



-S12 



-S23 



+B 2 Ls?j 1 ' l (s34, ti23\ S56, S12) + T j| m (si2, S34, S 56 ) 
+-B4 Ls^™ h (si2, *234; S56, S34) 



(13) (6 + |(l + 2)|3 + ) 
(12) [5 6] (4+|(2 + 3)|l+) 



(13) 



(6 + |(2+3)l|2-) Lo (^S") | [6 4](43) L °(=S) 



tl23 ^123 

, [4 2] (5~ I (3 + 4) |2-) 
J [43](65) (1-|(2 + 3)|4-) 

(5-|(2+3)4|3+) Lo (' f S) 



^234 



^234 



+ 



(2 3) 



(51) [12] 



^123 



[3 2] 



^234 



and 



^1 = -o 



(13) 



2 (12) (2 3) [5 6]*i2 3 (4+|(2 + 3)|l+) 



+ 



((6+|(2 + 3)12|3+))' 

+ ([6 4] (43) ti2s) 
1 [6 2] 2 



l 123 



(14) 



2 (12) [2 3] [3 4] [5 6] 
F s =F s i+flip(F s i). 

In these expressions, 

(13) ((6+|(l + 2)|3+)) 5 



B 5 = 



+ 



i>>2 



+ 



(12) (2 3)[5 6]ti23 (4+|(2 + 3)|l+) 

[12] 3 (4 5) 2 

[2 3] [13] (5 6)ti23 (1+|(2 + 3)|4+) ' 

(13) ((6+|(l + 2)|3+)) 2 
(12) (2 3) [5 6]ti23 (4+|(2 + 3)|l+) 
[12] 2 (45) 2 (2+|(l+3)|4+) 



(15) 



[2 3] (5 6) 1 123 (1+1(2 + 3) |4+) (3+ 1 (1 + 2) |4+) 
Bi = flip(B 5 ) , 

b 4 = mp(B 2 ) , 

and 

*123 (S56 + Sl2 — S34) — 2S12S56 



2ti 



B 2 

2 



(16) 



j [12] ((4-|(l + 2)(3 + 4)|5+)) -si2S 3 4 (4 5) 



2 (12) [3 4] (5 6) (1+|(2 + 3)|4+)(3+|(l + 2)|4+) 



T = Ti+fhp(Ti) . 

This amplitude is one of three leading-color he- 
licity amplitudes; in addition, five other primitive 
helicity amplitudes are needed to construct the 
subleading-color amplitudes for Z four jets. 
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